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Abstract 

The resonating mean-field theory (Res-MFT) has been applied and shown to effec- 
tively describe two-gap superconductivity (SC). Particularly at T = using a suitable 
chemical potential, the two-gap SC in MgB2 has been well described by the Res-Hartree- 
Bogoliubov theory (Res-HBT). The Res-HB ground state generated with HB wave func- 
tions almost exhausts the ground-state correlation energy in all the correlation regimes 
including an intermediate coupling regime. In this paper we make an attempt at a Res- 
MF theoretical description of thermal behavior of the two-gap SC. In the single-gap case 
we find a new formula leading to a higher T c than the T c of the usual HB formula. 
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1 Introduction 



A two-gap superconductivity (SC) of magnesium diboride MgB 2 with critical temperature 
T c = 39K has been discovered nearly a decade ago [JJ. Hitherto, intensive studies had been 
made to raise the T c of usual BCS superconductor in the weak coupling regime [21 [3j H] and 
to obtain Eliashberg's critical temperature in the strong coupling [3 El [7]. The T c = 39K is 
close to or even above the upper theoretical value predicted by the BCS theory [8j. Through 
ab initio density functional computations it has been estimated as 22K by Kortus et al. [9]. 
The existence of two energy gaps in MgB 2 has been predicted theoretically by Liu et al. 
[TO] employing the effective a and tt two-band model. They have obtained gaps at T = 0, 
A CT = 7.4 [meV] and A 7r = 2.4 [meV] and also their temperature dependencies and T C = 40K. 
The two-band model was first proposed by Suhl et al. [11] and next introduced by Kondo [12]. 

To approach such problems, Fukutome has developed the resonating Hartree-Fock theory 
(Res-HFT) [13] and Fukutome and one of the present authors (S.N.) have extended it directly 
to the resonating Hartree-Bogoliubov theory (Res-HBT) to include pair correlations [HI [15] . 
Steadily the Res-HBT has succeeded to describe effectively the two-gap SC [16] (referred to as 
I). If we get a Thermal Gap Equation in the Res-HBT, it is a strong manifestation of analogy 
of the Res-HBT with the usual BCS and HBT [2J El UU [TS] . The Res-HBT has the surprising 
fact that every HB eigenfunction in an HB resonating state has its own orbital-energy. Due 
to this fact, the Res-HBT also becomes a possible candidate and is considered to be useful for 
such subjects. This is because that the Res-HBT has the following characteristic feature: The 
Res-HBT is equivalent to the coupled Res-HB eigenvalue equations and the orbital concept 
is still surviving in the Res-HB approximation (Res-HBA) though the orbitals are resonating. 
This feature permits us to say that in some sense the band picture has a correspondence to 
the orbital concept in the Res-HBA though bands of different structures are resonating. The 
Res-HB ground state generated with HB wave functions (WFs) which are the coherent state 
representations (CS reps) [19], is expected to almost exhausts the ground-state correlation 
energy in all the correlation regimes including an intermediate coupling regime. 

To demonstrate the predominance of the Res-HBT for superconducting fermion systems 
with large quantum fluctuations over the usual BCS and Eliashberg theories, we already 
have applied it to a naive BCS Hamiltonian of singlet-pairing. A state with large quantum 
fluctuations is approximated by superposition of two HB WFs which are non-orthogonal CS 
reps with different correlation structures. We have optimized directly the Res-MF energy 
functionals by variations of the Res-MF ground-state energy with respect to the Res-MF 
parameters, i.e., energy-gaps. The Res-MF ground and excited states generated with the two 
HB WFs explain most of the magnitudes of two energy-gaps in MgB 2 . Both the large energy- 
gap and the small one have a significant physical meaning because electron systems, composed 
of condensed electron pairs, have now strong correlations among the fermions [161 120] - 

To go beyond the above ab initio density functional computations and phenomenologies, we 
develop a thermal Res-HBT. A thermal variation is made to satisfy diagonalization conditions 
for thermal Res-Fock-Bogoliubov (Res-FB) operators along a way very different from the usual 
thermal-BCS theory [211 [221 123] . We derive formulas for determining T c and the behavior of 
the gap near T=0 and T c . Particularly in the case of an equal magnitude of two gaps but with 
different phases, we find new formulas for T c boosting up T c to higher value than the usual 
HB's value and get new analytical expressions for gaps near T=0 and T c . In the intermediate 
temperature region, we solve numerically a thermal resonating gap equation. 
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2 Thermal resonating HB eigenvalue equation 



Consider the whole Res-HB subspace |^ Res(fc) ) = Er=i c^kXk); HB WFs), (Jfc = l 



n), 



in which the Res-state with index k : 



1 and the Res-states with indices k = 2,- ■ -,n stand for 
the Res-ground state and the Res-excited states, respectively. We introduce a temperature 
dependent Lagrangian with the Lagrange multiplier term E^ to secure the normalization 
condition (^ Res « | ) = l : 

^R^ELiE^i^^^ z=Z n r= i®^ (2.1) 

variation of which is made in a quite parallel manner to the one in the previous works where 

J>(r = l, • • • ,n) is a Res-FB operator [I3J E]. From the variation of with respect to 

(fe)* 



£™ =1 {tf [^Z] Res ]-£( fe )}-[det 



Cr "' for any k, we get a thermal Res-HB CI equation to determine ci 

Zrs]^cf ] = 0. 

We define D rs and D rs as D rs = u s z~l5u\ and D rs = 8u s z~lu\, respectively, where z r _ 
and uj = [bj,aj] (a r and 6 r :r-th HB amplitudes). The variations of thermal HB interstate 
density matrix H^[iT]R es and thermal overlap integral [detz] 1 / 2 are given by 



(2.2) 



8W rs \Z] Res 



S[det z r 



: D rs (l 2 N 

1 



■W rs [Z} Res ) + (l 2N -W rs [Z} Res )D r 



Tr(D rs + D rs ) ■ [det z r 



(2.3) 



We obtain also the variation of Hamiltonian matrix element i?[W rs [^]R e s] as 
5H[W rs [Z} Res ]^Tr{T[W rs [Z] Re8 ]5W rs [Z] Re8 } , 



J"[W rs [Z} Res 



rs;a/3 



[z\ 



Res 



F rs \Z\ Rcs 

-D* sr [Z} Res 
5H[W rs [Z} Rcs ] 



D rs [Z] R, es 



D 



rs;a/3 



\z\ 



Rcs 



Rcs 

_ 5H[W rs [Z] Res ] 
~~ SK* 



Kp + [a(3\^5]R rs . Sl [Z] Rcs , 



--[a-i\f35]K rs . 5l [Z] Res . 



(2.4) 



Following [H] and the Res-HF theory [13], writing L[Z]|| = X)Li E" s =i Cf s B{k) [Z} Res c r k> c s k) 
and Cfs^^Z]^ = {H[W rs [Z} Res }—E^}- [det Zrs} 1 / 2 , from the variation of L™ for any k, we 
obtain a thermal Res-HB equation to determine the thermal MF WF u r as 

ELi£:=iK-[iW^ fe) = o, 

K$$ [Z} Rcs ={(l 2W -H/ rs [Z] Res )^[H/ rs [Z] Res ]+^[H/ rs [Z]R CS ]- J E;( fe )} -W^W [det ^ r 
From (12.51) we can derive a thermal Res-HB eigenvalue equation which is expressed as follows: 

[Jv^ResWr]* = 6 ri U rh € ri = € rs - Y!k=l {H\W rr [Z\ Res } - E^} \c r k) \ 2 , 

F r [Z] Rcs = F r [Z} Rcs = F[W rr [Z] Res ] ELi I^P 1 (2.6) 



(2.5) 



+ELiE:= n i {ti$[z]^A k) *ci k) + JC^[z] Res c r k) ci k >} . ^ 



We call Jv[Z] Res the thermal Res-FB operator. For later convenience we denote W^s[i?]Res and 
•Tv^Res simply as W rs and T r . The Res-HB eigenvalue equation T r u rp = e rp u rp obtained in 
thus, can be extended to the thermal Res-HB eigenvalue equation (12. 6p in a formal way. 
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3 Thermal resonating gap equation 

The thermal HB density matrix is given as W„-p[.7>p] = [l+exp{/3J> p }]~ 1 (/3= 1/A;bT, r = l, 2) 
in momentum p. Using a Bogoliubov transformation gi(2) p , W / ii(22)p[^ r i(2) P ] are diagonalized as 



rp 







1 -w. 



rp 



l + e@ €r P 1+e pCr P 

which are the generalizations of the usual results to the Res-MFT. This is proved elsewhere. 
By making the Bogoliubov transformation g rp , eigenvalues ? rp are obtained by diagonalization 
of the thermal Res-FB operators T rp with additional terms (H[W rrp \— E)\c rp \ 2 (r = l,2). The 
thermal HB interstate density matrix in the whole Res-HB subspace is given as the direct sum: 

W P [F P ] =g P W p gl = @ 2 r=1 W rrp [F rp ], W rr [F T ] =g rp W rp gt p . (3.2) 

Suppose a tilde thermal Res-HB density operator Wu2) p for equal-gaps to be 



l — w. 



rp ' 



,-/3 e r 



(3.1) 



W 1(2)p : 



W{ 







I T 

mi 12 



tu) . 

1(2)P ' 



W l(2)p J 2 







(1-4%) - /2 



(3.3) 



Here J 2 is the two-dimensional unit matrix. Performing the unitary transformation by <?i( 2 ) P ! 
we obtain the following thermal Res-HB density matrix W^) p - 





K1/ l(2)p - »l(2)p ^ K l(2)p#l(2)p - yi(2)p 



W l(2)p J 2 







:i-4 ( 2)p> j 2 



a W)t 

yi(2)p 



(3.4) 



i (1-cos 1(2)p (l-2^y } I 2 {-(+)} x |sin e 1(a)pC -^w (1-2^ ( %J • h 

_{-(+)} x Isin 1{2)p e*( 2 ) (l-2wl^ p ).I 2 \ {l^os 1(2)p (l-St^^)} h 
In the equal-gaps case (iffWn] = i/|W 22 ] = -^[H / ]), following I, the Res-HB ground (excited) 

pRes 
v gr(ex) 

Case I 



energy E^Jf is classified into two cases, according to the solutions for the Res-HB CI equation: 



E R f , 

gr(cx) 



Case II 
E R f « 

gr(ex) 



#[W/]-if[W/ 12 ] > 0, 

' -•fe[H/]+(-)F[H/ 12 ]-[det^ 12 ]5 



l+(-)[detz 12 ] 



#[W]-F[W 12 ] < 0, 

i -r-fe[WH- (+)H[W 12 \- [det^ 12 ]l) . 

1- (+)[det z i2 j2 v 7 



(3.5) 



(3.6) 



The Res-FB operator for spin-up state, accompanying quantities with upper or lower 



sign corresponding to Case I (13.51) and Case II ( 13.61) . is expressed as 



(2)P" 



Ke^h {+(")} ^I,// 2 

{+(-)} xTjp.Ja 



rt 
Ap 



(3.7) 



The quantities JijJ ej] , and for Case I(upper sign) and Case II(lower sign) are defined as 



sin 2 |(cos 2 f 
cos 9 n 



A 2 

=F [det2;i 2 ] 

En 



l±[det^i 2 



Tl El es = H[W}-El es ) (3.8) 



p - 



where 



-Tl=-h\^N(0)V-arcsinh ± [detz 12 ]tj 



A 



[det z 12 ] 



exp 



D. 



2N(0)hw D { ln(l + x 2 )+2x-arctan| 



(3.9) 
(3.10) 



At finite temperature, the quantities J+ e , 7"_ Ej ,, J 7 ^ and A become temperature-dependent. 
This is explicitly expressed by a subscript T. Using the distributions (13. ip we require corre- 
spondence relations cos 6 pT =>■ cos 6\ p and sin 6 pT =>• sin 6 lp given through 



cos 6 r< 



4- T 

' r +e v T + • r - 



sin 9 T 



e 2 +Af 
A T 



2 e\ p 



A 



'7 1 



(l-2«$ p ) 



-255^ 



• r +e p T + ^-i 



(3.11) 



■2w\ p 



Notice the multiplication factor l—2w^ p \ The e£ p (=e|) is the quasi-particle (QP) energy: 

^={(^ £p T+^-£ t)^ • The first two equations in parenthesis of (13. lip are unified into 
the second single equation in (13. lip . From now it is shown that equation (13. lip plays the role 
of the self-consistency condition at finite temperature. Dividing numerator and denominator, 
respectively by (e 2 + A|,) 3//2 in R.H.S. of the second equation in (13. lip , we have 



1T\ 



( £ 2 + A 2 )I \ A T 



(1-255+) 



S 7 



rt 



1-255! 



w 



■ w 



lp- 



(4+ A f) f 

Now we demand a new condition for Thermal Gap Equation 



1 f jrt , jrt 

- [ - r +e p T ^ —£ P T 



(4+A 2 )t 




(l-2«5) = 0, 



(3.12) 



(3.13) 



which leads to 

l-iy(0)V- arcsinh 



Cr) 



+ E^ {±) hw D j: p B p TA 2 T ■ [det z 12 ]f £ p C p = 0, 



Res(±) 



gl'T 



D 



■ ±2/V(0) hiODXrp- arcsinh 



inh ( — J 



2 -JV(0)V- arcsinh 



[det z 12 ] 



± [det z 12 ] 



(3.14) 



which we have calculated using the solutions for the Res-HB CI equation obtained in I. We 
also give the following definitions for ^2 p A p , J2 p B P and J2 p C p : 



E P A" Ep-^p' 2~2 p C P — Y2 P 



(1-2551 



(3.15) 



L (4+ A T) f 4+ A t (e 2 +A 2 )f_ 
Rearranging (I3.14p . it is cast to a Thermal Gap Equation with similar form to the one in I: 

i 



N(0)V 



arcsinh 



l±2iV(0)- 



arcsinhl-j A r -A r £ H 



1 ± [det z±2 



-•[det Zi 2 ]| 



l+<TlT 



:4/V(0)- 



arcsinhf- 



{^) Ar-A T Z P B P 



■ [det z 12 ] 



(3.16) 



l±[det z 12 ] 

The summations J2 p A p , X]p-^p anc ^SpCp near T = are computed in Appendix A. Substituting 
( 1A~TT) . ( IA~T2l) and ( lATH into ( l3~T6|) . then, near T = we have the gaps for Case I ([33]) as 
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1 - 



[ det *i2]o J arcsinh f^A _ i , [detfu]! \ . T i 
N(0)V arCSmh UoJ 1+ N{0)V T \ 



(3.17) 



On the other hand, with the aid of [det 212] o — 1 — 27rA r (0)/ic<;£)a;o(0 < xo 1) which is easily 
derived from Taylor expansion of (I3.10p . we obtain the gap for Case II (13.61) as 



1 + ; 1 ( 7T - (l-N(0) V-y/l + (N(0)V) 2 ) arcsinh ( ^) } ■ T? 



7Ty/l + (N(0)V) 

where we have used the approximate relation 



,(3.18) 



— N(0) Uarcsinh jarcsinh — 1 — lj* ~ arcsinh f /±J 



-1. (3.19) 

In ( I3.17P and (I3.18p . as shown from Appendix A, the quantity rf 11 ^ is defined by 

r ? «> s ^ {Ar}" ? {g}'«P (-0Af») .(-M,..) (3,0) 

In the opposite limit T— (T c for Case I) the gap becomes very small, [det z^]^ 2 — >T, 
then J r 2 T ^-AiV(0)Varcsinh(^; i? /A)/4 and (J 7 ^ T +.F_! e T )/2^e p /4 if we use the relations 

(1378]) and ([379]). We have an approximate QP energy e^ 1 ~ [^+{-AA^(0)Varcsinh(na; D /A)} 2 ] 1/2 /4 
in which the appearance of numerical factor 1/4 should be become aware carefully. This is 
because two HB WFs have different correlation structures ip2 = 7r and ipi = 0. In such a case, 
returning to the original form of the BCS gap equation but with the modified QP energy e^ 1 , 
the thermal gap equation is expressed as l=V/2 J2 p {^~^p) l^p and leads to the integral form 

U^lVjTy^y. (3.21) 

We should emphasize that this form results from the above numerical factor 1/4. Intro- 
duce a dimensionless variable y],^e/8k B T l , its upper-value y, 1 l =hjj B >/'8k B Tl and the Debye 

c 

temperature Q B) = ftu3 B) lk B . Integrating R.H.S. of ( I3.2ip by parts, it is approximated as follows: 

1 , t r°, , l2 1 t , f^ c \ , fe c hw D \ , (e D \ 

— — ~ln?/ir-/ ay iny seen y — layL,r+m\ =ln — - — — - = in — , (3.22) 

4iy(0)V Jo T * V * J Vitk B T l J \Tl) 

where we have used the formula given in textbook [24J. The number C is the Euler's constant 

(C~ 0.5772) and e c ~ 1.781. Finally a small rearrangement yields 

T* = O.2830£>e~*^. (3.23) 
It should be compared with the Eliashberg's formula [23] and the usual HB's one for T c [S] 

T c = 1.130# D e~W\ (3.24) 
The new formula (I3.23|) gives a high critical-temperature, e.g., T C : = 72.87K for iV(0)V=0.25 
and #d=700K. This T* is in contrast to T c obtained by the famous HB formula (I3.24p given 
by Rickayzen and Cohen in (8], namely, T C =14.49K for the same values of N(0)V and 9d- 

From now we discuss behaviour of the gap near T c . In the above the modified QP energy 
e=[e 2 +{— AA^"(0)l / arcsinh(/ici;£)/A)} 2 ] 1 / 2 /4 plays a crucial role to boost the T c ( I3.23|) comparing 



with T C = 14.49K by (I3.24p . Note the numerical factor 1/4 in e. First consider A^ near T^. 

Using this form of the modified QP energy, the gap equation is roughly rewritten as 

1 f huJD , 1 1 ( e \ 1 ( £ \ 

- I cfc-tanh — — — ~ / cfe-tanh — — — 
4 Jo e \2k B T) Jo e \8k B T ) 



AN(0)V 4 

- ( ArN{0)V- arcsinh f ^ ] \ [ °de { \ tanh ( -i— ) - ^--±— sec h 2 



J) Jo U 3 \8k B TJ e 2 8k B T \8k B T 



(3.25) 
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from which we obtain 
1 



iV(0) Vxt • arcsinh 




2 



N(0)Vxt • arcsinh 



(3.26) 



where hwo/ksT = e c /2tc ■ hwD/ksT and we have used the famous integral-formula (1A.18I) 
Using arcsinh (1 /a; T ) ~ In (2/x T ) (0<x T < 1), (E22D and (EEZBjl . we get A^ near T\ as 



AL~2ttJ , , f ° 1- — j W c r • (3-27) 
V 7C(3) N(0)V \ Tl )\j Tl K ] 

Such a formula has been brought through the use of the modified QP energy e which owes 
to a resonant feature of the multi-band SC. This new formula shows a more complicated 
temperature-dependence of than the a/^T a/TJ — T dependence of A^ presented by equation 
(36.6) in textbook |22j and by equations (16.32) and (16.33) in textbook [23], respectively. 
Next, for Case II, (J^ £ T + J^ £ T ), j£ and become infinite simultaneously in the limit 

1 /2 

Ay— >0 (xr— >Q) due to the existence of 1— [det zu}^ in the denominator. Then mathematical 
handling for such a problem is too difficult and therefore we can not easily get a formula for 
T c n in an analytical way as we did in Case I. Let us denote T c for Case II as T c n . At T~T c n , 
Aj> almost vanishes and 1 — [det Zy^j 2 — >2ttN(0)HujdXt- Using the J-\ expressed as (I3.9p . 
we reach to the following asymptotic forms: — >•— (47rA ir (0)) _1 A ir (0)V A arcsinh (Hud/At) and 
( 4 £pT +P £pT )/2^(A7rN(0))- 1 e P /A T . The QP energy {(4 £pT +P £pT ) 2 /A+^ 2 T }^) is approx- 
imately calculated as ?J II ~(47rA ir (0)) _1 £:p/A T (0< A r ^C 1). Here we discard the contribution 
from comparing with the contribution from (J* T + P £ T )/2. Returning again to the 
original form of the BCS gap equation but with another modified QP energy e£ n , the ther- 
mal gap equation is obtained as l = V A /2^ p (l — 2w^)/e£ n which also leads to the integral form 

1 Wo £ V2A: B T II -47rAr(0)A^V 

We introduce the dimensionless variable t/^ =e/(2/c^r-47rAr(0)A^) and its upper- value y^u = 
hu D /^k^T u -4:7TN(0)A^). Integrating fl3~28l) by parts {e D /k^T ll = e c /27i-hu D /k^r 11 ), we have 

4^^V^" ln * + ln (^) =ln (^W^) ' ^'^ 

which reads 

&d | 1 I #D 



At 7riV(0)T n ^ i 4iV(0)y7riV(0)A? j ~ ^#(0)7* V 4iV(0)y7riV(0)A? ) ' (3 ' 30) 
From ( 13 .301) we obtain a quadratic equation for A^ very near T c and then we have a solution 

A 11 = 1 6 ° 1 (* Q1) 

T 7rAT(0)f n 4AT(0)V7rAr(0)' 1 J 

in which at T 11 = T c n , the A^ vanishes. Finally we can determine the critical temperature Tj 1 

for Case II as c 

Tl 1 = —0 D N(0)V = 1.3340 D iV(O)V. (3.32) 

71 



The simple formula (13. 32ft gives a high critical temperature, e.g., T c n = 198K for N(0)V = 0.25 
and 6d = 700K. Finally A^J near T c n can be approximately obtained as 

A "~- e -^^^l (3 33) 

which is linearly dependent on T-T^ 1 . It is very interesting that we could find such a dependence 
of A 1 ^, comparing with the usual dependence a/ T— T c n of A 1 ^. The numerical results for N(0)V= 
0.25 and 6fc>=700K, obtained from 1ET25]) . (l3TMjl and (KT52]I . are illustrated in Figffl below: 

T Critical Temperature T c 

1 98K • 



72. 87K 



14. 49K 



Figure 1: Critical Temperature T c : 1. T C I = 72.87K; 2. T C = 14.49K; 3. T c n = 198K 
In the intermediate temperature region, substituting (lA.16j) and (IA.17j) into (I3.16p . we have 

c 1 \\{ , . . fl\V 2tt 2 (e c 1 

m 



^ ±[d e t2l i 4 f (0)yarrail Kiyj 2K(3) 



7 D 



X 



jiV(O) Varcsinh (J-) -(l- [det 2 12 ]| 



2 /e 



7T 



1 



l±[det z i2iT 

2 



TV 



J (3.34) 



3 \ T l±[detz 12 ]|y 



x|.[det 2i 2 ]£. 



A solution for equation (I3.34p is classified into the following two cases: 

Case LThe R.H.S. of (I3.34p is approximated to be zero since < xt <C 1, from which we 
obtain an equation to determine At for a given T as 



N(0)V 



— Xyarcsinh I ■ 



(1 



2tt 2 



n l+[det* 12 ]f W V 21 ^ 3 )\ 

Using the relation and the approximations 



■jo 



In I 



1 



9 



D 
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l + [det zi 2 ] 



T 



T 



arcsinh ( — 

\X 



x V x I x V ar x 2r 



(3.35) 



(3.36) 



and e x ~l+x + ^x 2 , for [detz 12 ]^ ~0.3 we have 



1 

XT + -- 



V0J82 







T\ s 



from which, finally we have a very simple solution for x*r (A/"(0)V = 0.25 and 6 I £> = 700K) as 

(3.38) 



y/0~7&2 / , . f T\f T 
'ln(0.436)-ln 



0.436x0.25 



700/ V70Q 



8 



1/2 

Case ILUsing 1— [det 2i 2 ] T m2irN(0)huDXT, from ( 13.34(1 we have 



(0.283) 2 {iV(0)V} arcsinh 2 (-^)-0.782{7riV(0)^; D } ; 



ln[ 



0.283 



— ln( — ) + arcsinhf— 



T 



T 



(3.39) 



27rN(0)hu D J 

|iV(O)Varcsinh0^ - (l-[det z l2 ] T H = -^(0.283) 2 [det z 12 ]\. 

Expanding (I3.39P with respect to arcsinh(2/xT) and neglecting a constant term which is very small 
for [detzi 2 ]T /2 ~0.3 and for (/V(0)l/= 0.25, N(0)hu) D = 0.01 and Q D = 700K), (I3T39]) becomes to be 
a quadratic equation for arcsinh(2/a;x). Finally we have the following solution for xt'- 



xt = 2/ sinh 



1 Jl— 0.3 0.782{ttx0.01} 2 /T x21 

2 I 0.25 + (0.283) 2 (0. 25) 2 V700/ 



+ 



\ 




0.782{vrx0.01} 2 



0.782{vrx0.01} 2 



(0.283) 2 (0.25) 2 \700/ J (0.283) 2 (0.25) 3 [ 



0.25<Hn 



0.283 

27TX0.01 



■In 



T 



70Q 



TV 
700/ 



We draw below the numerical results of the solutions for Case I and Case II. 



(3.40) 



x t Temperature Dependence of the gap x T for Case I 
0.25 

0.2 



Temperature Dependence of the gap xr for Case II 




250 




Figure 2: Temperature dependence of the Figure 3: Temperature dependence of the 

gap, Case I for [det z l2 } 1 ^ 2 = 0.3, N(0)V = gap, Case II for [det ,z 12 ] T /2 = 0.3, N(0)V = 

0.25 and f5 Z) = 700K. 0.25, N(0)hw = 0M and fl D = 700K. 

The formula for Case I gives a high T c I ,e.g.,T c I =72.87K for parameters iV(0)^=0.25 and 
#£>=700K. This is in contrast with T c of the usual HB formula giving T C =14.49K for the same 
values, A^(0)y = 0.25 and #d = 700. The formula for Case II gives also a very high T c n , e.g., 
T ( J I =198K for the same values of the parameters. They are illustrated together in FigJTJ The 
temperature dependence of gap near T=0 and T c becomes more complicated than that of the 
HB and Abrikosov's descriptions j22j [23] . At intermediate temperature, as shown in Figs. [2] 
and El we have got the solutions of At for Cases I and II. We assume [detzi 2 ]r «0.3 (Cases I 
and II) and N(0)hu) D =0.01 (Case II) to acquire real solutions. Anyway we could obtain really 
the solutions xt=0.030 (Case I) and 0.243 (Case II) for T=80K. The former has a negative 
gap below 72K. To our great interest, that value is almost equal to the T c given by (13.231) . It, 
however, recovers a positive and small gap. Further it increases as temperature rises up to 
around 190K but shows vividly a decreasing tendency beyond around 200K. In this sense the 
former is considerably good solution. On the contrary, the latter solution naturally decreases 
to 0.198 as temperature rises up to around 2000K but never vanishes. This means the latter 
solution has no tendency approaching the T c given by (I3.32p . Much improvement of the above 
results should be possible if the original equation ( 13.34(1 can be solved more accurately. 
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4 Summary and further perspectives 



In the present paper we have made an attempt at a Res-MF theoretical description of 
the thermal behavior of the two-gap SC. To show the predominance of the Res-HBT for 
superconducting fermion systems with large quantum fluctuations over the usual BCS and 
Eliashberg's theories, we have applied the Res-HBT to the naive BCS Hamiltonian of singlet- 
pairing. We have derived gap equations within the framework of Res-HBA. From the Res-FB 
operators T\ and Ti with equal-gaps, we have found the diagonalization condition for them, 
which is essentially the same form as that of the BCS theory. It leads to the self-consistent 
Res-HB gap equation, from which we could derive the present gap. Here we have concentrated 
on the derivation of the Thermal Gap Equation with the use of the thermal Res-HBA. From the 
thermal Res-FB operators Trr and J~2T with equal-gaps we also have found the diagonalization 
condition, which is just the same form of the condition used at T = 0. This reads the self- 
consistent Res-HB Thermal Gap Equation and makes possible the derivation of the present 
new formulas to determine the T c and the gaps near T = and T c . 

For unequal two-gaps, it is also possible to realize the above diagonalization condition 
for Res-FB operators J- rp (r = 1,2). Transforming by a unitary matrix g rp , T rp is easily 
diagonalized. Noticing the same correspondence as the correspondence in (13. lip . cos# rp =>- 
cos6* rp and sin rp =^sin 9 rp , we assume each diagonalization condition ( 13. 12ft holds even in this 
case. Then we obtain coupled equations through a function of Ait and A 2 t expressed as 



f 2 / 2 F f 



{el + A 2 r )3/2 I A r 



rT 



(i - say 



(e 2 + A 2 T ) 3 / 2 \ rep rep ) \ ZWr Pj 

which reduces to equation in R.H.S. of ( 13.1 ip if A 1t = A 2 t- The quantities J~l& T and ^t± £p 
are given by the equations similar to (5.9) in I but with more complicated forms of A^ and 
A.2T- For the time being, as was done in the previous section we here also use the function 
(£p + A 2 T ) 3//2 by which we divide numerator and denominator, respectively, in (14. ip . After 
equating the numerator to the denominator and using the relation 1 — 2w^ p = tanh(e r /2/cfiT), 
we sum up over p, namely integrate both sides of the equation over e, to achieve the optimized 
conditions. Thus we obtain Res-HB coupled Thermal Gap Equations and reach our ultimate 
goal of computing thermal two-gaps. Along such a strategy and method, we will make a 
numerical analysis to demonstrate the thermal behavior of two-gaps. 

To describe a superconducting fermion system with N single-particle states and to approach 
such the problems, we must provide a rigorous thermal Res-HBA and MF approximation. We 
have an expression for a partition function as Tr(e~^ H )=2 iV ~ 1 f (g\e~P H \g)dg in an SO(2N) CS 
rep \g) [TJJ] where fdg is the group integration on group SOfiN). Following Fukutome [TB] . 
introducing the projection operator P onto the Res-HB subspace, the partition function in 
the Res-HB subspace is computed as Tr(Pe _/3H ). This is a purely theoretical object and can 
be calculated within the Res-HB subspace, e.g, by using the Laplace transform of e~^ H and 
the projection method. The result leads to an infinite matrix continued fraction, a concrete 
computation for which, however, is very difficult. On the other hand, a calculation of the 
Res-HB free energy is carried out in the same way as usual. This is made parallel to the usual 
thermal BCS theory, which will be given elsewhere in a separate paper in a near future. 
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Appendix 

A Calculations of ^2 p A p , ^2 p B p and ^2 p C p at zero and in- 
termediate temperature 



First, equation ( 13. 16ft is shown to reduce to the Res-HB gap equation (4.10) in I as T— >0. 
Using a variable e = £,At instead of s, summations ^2 p A p , J2 p B p and ^2 p C p near T = are 
computed as follows: 



%M = arcsinh (-) ~ -7n==^f + A ^ A ^ = - T ' (n) + • ' ' > ' 

2JV(0) \x T J ^Jl + x% 2 

At E p B p 



2N{0) 
2JV(0) 



arctan ( — ) + B(T), B{T) 



1 2 I 



,1(11) 



1 



+ c(t), c(t) = -t; 



2 

I(II) 



T I(II) 



(A.l) 



detailed calculation of which is given below. With the use of the relations (I3.9P and Ar/e p - 
— 2J^ T / (J] eT +J^ t ) which lead to e = a/£ 2 + 1 /wuAy , E P A> i n ( 13. 15|) is converted to 

E P A> /"* „ £ 2 n f x £ 2 1 .,iuii_--m-.i 2) 



2JV(0) 



u (e+i) 



l+e l 



,1(11) 



where d I( W ^> 1 for large frhJo and for Case I and Case II, A^ is defined as 



A 



i(ii) . 



, = -x T < N(0)V- arcsinh — + (-)[detzi 2 l 
' 2 T \ y } \x T ) 1 J[ 12J J l+(-)[det^ 12 ] 



1 



(A.3) 



Introducing a new variable y = a/£ 2 + 1, (1A.2|) is integrated partly and approximated to be 

dy ; e d ^ + 2/ dy— , e~ d ^. A.4 



P - ~ arcsinh f — | - 



2JV(0) 



v9 



Similarly, we get approximate formulas for ^ B p and E„ C p as 

a, v 



'P / 1 \ f 1 1 1 

— — ^- — ~ arctan — — 2 / ct-u == 

2iV(0) \x T ) A 



1 



y 



(A.5) 



(A.6) 



2iv(o) yrT4 k y 2 Vf^i 

To carry out integral calculations in (1A.4j) ~ (1A.6j) . it is convenient to use an integral repre- 
sentation of Bessel function The Bessel function of order v is represented as 

/£(§)* 



dy(y 2 - \y-2e 



Then, the integral form of the Bessel function of order and its exact result are given by 

^ f°° . l „. A¥ 



(A.7) 



Kn(z) 



r(§U dv 'Vy^ 



-zy 




2z 



e~\ (r (I) = VtF)- 



(A. 



Using (lA.8p . integral calculations of ( 1A.4j) ~ ( 1A.6j) are made in the following ways: 



dy . 

i y^F^i 



-dy 



vy 



i 



dze~ zy 



1 Jd 



dzK (z) 



,_I 1 ,_1 3 5 
d 2 d 2-| d 2 

2 4 



(A.9) 
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i; 



dy 



-dy 



y 2 Vi 



dy 



dwe- wy 



1 Jd 

, 1 __3 9__ 

d 2 - d 2 + -d 

4 



d 
d 



dw I dzK (z) 



As a result, we obtain the approximation for A(T), B(T) and C(T) near T = as 

A(T) ~ 0, B{T) = C{T) ~ -lf n) . 

2 

Further in the above near T = we make the following approximations: 

1 



arctan 
arcsinh 



7T 

-2 -x T , 



1 - £ £t, (0<io<l) 



/1\ 1 ,- - 

arcsinh I — I (xt — Xo), [det z^jy — [det -2i 2 ]^ =0 . 

V''o/ •''(! 

1/2 — 

Next, let us introduce a new variable y by e=4(lzt[detzi2] T )ksTy and quantities xt=^t/^d 
and y T ±) = y^ 2 + A|{4(l ± [det z l2 ^ 2 )y l hjj D lk B T where A T = A T A(0)\/arcsinh (hw D /A T ). 
In intermediate temperature region the modified QP energy e is approximated as = 

£ 2 +A 2 ,{2(l±[det z^r 2 )} -1 - ^ £>A T , EjA> E p 5 p an d E p C P in are recast to the 

following integrals up to A^: 

A t>huj D -I 

2.„ A> _ /• & _ 



(A. 10) 



(A.ll) 



(A.12) 



2A(0) 



tanh 



" ,/ c 2 



|(±) 

21cbT. 



A 2 T 



T 



de- 

e 2 \/e 2 +A 2 



tanh 



|t±) 
2k bT. 



dyil 



y { T ] x T 



=tanh 



y 2 +[y { T ] XT 



(A.13) 



A(0) 



de- 



=tanh 



-hu D £ \ e 2 +A 2 



i r** 13 , A3 

. t anh 



2k B T) 2j_ hWD 3 / 2j A2 



eW+A 2 



g±) 



(A.14) 



2A(0) 



~ ' de- 

'" £ 2 \/e 2 +A 2 



tanh 



2fc B T, 



^ A 4 
& tanh 



£ 4 A /£ 2 +A 2 



£(±) 

2fc B T, 



(^t) 2 ^ J 1 - 



y { T ] x T 



v'W y z +[y { T ] XT 



; tanh 



y 2 +[y { T ] x T 



Further equations ( 1A.13I) and ( 1A.15I) are approximately computed, respectively, as 



2N(0) 



de- 



e 2 +A T 



~~rr tanh 



2A(0) 



de- 



-rrtanh 



e 2 +At 



71 



21 C(3)fe ±} JT 



2tt 



7 



,2A;rT/ tt 



C(3)(|/^ r ) 2 



(A.15) 



(A.16) 
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Taking only a leading term, finally terms A p and C p in ( 1A.16I) are approximated to be 



2N(0) 



-,C 



In 



17 1 ± [det z n \l T 



a 



±) ^VT^pCp 2 



(±) 

T ' 



2JV(0) 



Irp , 



„(±) 



21C(3) fe c 



± [det z 12 ) T/ 
To derive (IA. 16[) we give a integral formula 



2 / ~ \ 2 

k B T , 



(A.17) 



dy { — tanh y -seen 2 ?/ 

2T 2T 



7T 



rC(3), 



(A.18) 



which can be derived by using the famous mathematical formulas 
1 - _™ 1 

y 



tanhy-SE^x^^^^, y3 



- tanh y - ^sech 2 y = 64^^U 



{{2m-l) 2 7r 2 Uy 2 Y 



(A.19) 



Adopting a new integral variable y— (2m— l)7r/2-tan#, an integral of the second formula in 
( 1A.19j) is easily carried out for Hud ^> 1 as 



64 



oo 

7TS=1 



: 32£: 



1 



d0- 



7 



{(2m-l)% 2 +4y 2 } 2 ^ m=1 (2m-l) 3 7T s J 1+tan 2 tt 2 
where we have used E^iO™- 1 ) -3 = ( 7 /8) <(3), and C(3) =tt 3 /25. 79436 



C(3), (A.20) 



To get a finite value of ^2 p A p , expanding (1A.13I) around y^XT, (1A. 13h is boldly approxi- 
mated as 

- / dy- tanhy (y^xj) [ dy l\ tanh?/ - \sech 2 y\. ( A -21) 
2iV(0) Jo y 2 \ J J [y d y 2 J 



In a similar way we also get a roughly approximated integral form for ( 1A.15|) as 



2 



{ — tanh ?/ -sech 2 ?/ 

y 6 y 



(A.22) 



Integrations of (1A.21j) and (1A.22I) are easily made by using the integration formula ( 1A.18j) if 
we take the upper-value y^ to be infinite. 
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